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Abstract. Using atomic force microscopy on silica and float glass surfaces, we give evidence that the 
roughness of melted glass surfaces can be quantitatively accounted for by frozen capillary waves. In this 
framework the height spatial correlations are shown to obey a logarithmic scaling law; the identification 
of this behaviour allows to estimate the ratio kTF/iry where k is the Boltzmann constant, 7 the interface 
tension and Tf the temperature corresponding to the "freezing" of the capillary waves. Variations of 
interface tension and (to a lesser extent) temperatures of annealing treatments are shown to be directly 
measurable from a statistical analysis of the roughness spectrum of the glass surfaces. 



PACS. 68.35.Ct Interface structure and roughness 
related phenomena - 68.37.Ps AFM 



1 Introduction 

Glass is becoming a widespread substrate in high-tech 
developments like flat displays ,1, thin foiled x-ray tele- 
scope mirrors and lithography masks |2J . In all of these 
applications, controlling the glass surface roughness is a 
critical issue to enhance performance. Glass manufactures 
spend considerable effort in optimizing their forming pro- 
cess. The introduction in the late fifties of the float pro- 
cess to produce flat glass represented a revolution in the 
glass industry: the glass sheets obtained in this way are 
almost perfectly flat so that the expensive polishing step 
revealed to be unnecessary for most applications. The pro- 
cess somehow mimics and extends the traditional method 
of fire polishing. The latter consists of first heating the sur- 
face to a temperature high enough to obtain local melting 
and then quenching the resulting smooth surface. Laser 
polishing 0] can be seen as another modern extension of 
the same strategy. In the float process, the molten glass 
flows from the furnace on a liquid tin bath. The (liquid) 
glass thus spreads out on the tin surface to form a ribbon. 
The progressive cooling imposed along the tin bath allows 
the two glass surfaces (tin/glass and glass/atmosphere) to 
freeze without external mechanical contact. Note that sim- 
ilar comments could be made about other glass forming 
processes like fiber drawing. 

The residual roughness of these "fire polished" surfaces 
is extremely low, typically below the nanometer range 
when measured over a micrometric scale and has been 
attributed to the freezing of the capillary waves of the liq- 
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uid interface in the glass transition temperature range[S]. 
This particular phenomenon has recently been studied on 
glycerol HJ|7| and supercooled polymers[Sj by X-ray reflec- 
tivity. Though the surface roughness induced by capillary 
waves is long range correlated, the X-ray reflectivity of 
a liquid surface was shown to depend mainly on the to- 
tal interface width a oc fcT/7 2J, where k is the Boltz- 
mann constant, T the temperature and 7 the interface 
tension. When varying temperature in the glass transition 
regime, a clear saturation of the interface roughness was 
obtained below a threshold temperature |S][7||S]. A satis- 
factory quantitative agreement with a capillary waves de- 
scription is obtained in the case of polymers; for glycerol, 
however, a density driven layering effect was suspected, 
only a qualitative agreement could be obtained. 

In this paper, we present a detailed quantitative atomic 
force microscopic (AFM) study of the roughness spectrum 
of melted glass surfaces, namely silica surfaces obtained af- 
ter various annealing treatments and bare industrial float 
glass surfaces. AFM has previously been used to study 
the roughness of glass surfaces 10;; but only recently, a 
few reports (1 11112) have been given of a qualitative agree- 
ment between the values of the roughness measured on 
glass surfaces and the simple estimates obtained from a 
description of the surface in terms of frozen capillary waves 
a w kTc/j. Here again, a is the standard deviation of the 
height fluctuations along the surface (RMS roughness), 
7 the surface tension of the glass and the glass transi- 
tion temperature Tq is used as an estimate of the freezing 
temperature of capillary waves. Focusing on two differ- 
ent glass surfaces (silica and float glass) we show here 
that the height correlations along the surface give a di- 
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rect and quantitative access to the frozen capillary wave 
spatial spectrum. The use of these two different materials 
allows us study two different aspects of the phenomenon. 
Performing heat treatments on silica at different tempera- 
tures in the glass transition range we can study the effect 
of thermal history on the glass surface roughness. Per- 
forming measurements on the two faces (namely "tin" and 
"atmosphere") of a float glass, we can test the effect of the 
two different interface tensions due to the asymmetry of 
the float process. 

The paper is organised as follows. We first recall some 
results on capillary waves at a liquid interface and we give 
quantitative expressions expected for the height correla- 
tions in this framework. Then we detail the experimental 
protocol and we finally present the AFM measurement re- 
sults obtained respectively on silica and float glass. Finally 
we discuss their interpretation in terms of frozen capillary 
waves. 



2 Capillary waves 

Assuming a constant value of the surface tension 7, the ex- 
cess of energy E(h) due to the height fluctuations h(x\ , x%) 
of a liquid interface can be estimated to: 



E(h) = 7 / (Vl+lV/i 2 !-! d Xl dx 2 



\Vh\ 2 dxidx2 



(1) 
(2) 



where x\,X2 are the horizontal coordinates and 7 the in- 
terface tension. Using a Boltzman distribution, we easily 
obtain the expressions of the RMS roughness a and of the 
spatial correlation function G(r) and of the height vari- 
ogram H(r): 
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where H? = xiel + x-i~&\ and r — \~J*\. As the expres- 
sions are logarithmically divergent, an upper (Aa/) and a 
lower (Am) spatial cut-off must be introduced. Gravity ef- 
fects are neglected in the above description and the upper 
cut-off is given by the capillary length l c — y^f/pg where 
p is the glass density and g is the gravity constant. For sil- 
ica (7 ~ 0.3J.m~ 2 , p ~ 2.2kg.m -3 ) this length is approx- 
imatively equal to 4 mm. The lower cut-off is estimated 
to a molecular length £ . For silica and silicate glasses in 
general the structural unit is the SiC>4 tetrahedron, and 
we can take £q — 0.5nm. The typical roughness of a sil- 
ica surface in the glass transition region (Tq ~ 1450) due 
to capillary waves can thus be estimated to er w 1.5nm. 
In a first approximation we consider that the interface 



height fluctuations freeze at this glass transition tempera- 
ture Tq so that equations ll'ill can be used to describe 
the roughness of the glass surfaces at room temperature 
simply replacing T by Tq. A detailed treatment of the 
freezing of capillary waves based on the visco-elasticity 
theory can be found in [0]. 

Molecular and capillary length thus define the ulti- 
mate physical bounds of the expected logarithmic scaling. 
However, in the framework of roughness measurement, the 
cut-off lengths are imposed by the limited bandwidth of 
the experimental set-up. In the context of atomic force mi- 
croscopy the scan size, usually in the micromemeter range, 
defines the upper cut-off \m and the nanometric tip radius 
(or the sampling length) define the lower cut-off A m . 



3 Experimental method 

We detail here the sample preparation (annealing treat- 
ment for silica and cleaning procedure) and on the rough- 
ness data extraction from the raw AFM data. 



Sample preparation 

Float glass 

Float glass is a standard industrial product with approxi- 
mate composition: Si0 2 72%, Na 2 14%, CaO 9%, MgO 
4%. Note that due to the asymmetry of the float process 
some traces of Sn are present in the vicinity of the "tin" 
surface. 1 cm 2 samples were cut from a 3 mm thick glass 
sheet. No additional treatment but cleaning (described be- 
low) was performed on these surfaces. 

Silica 

The silica samples are commercial silica wafers (vitreosil, 
EQ512), 2 mm-thick, two side polished. Square samples 
of 1cm x lcm were prepared. Prior to the annealing treat- 
ment, the sample preparation was as follows: the sample is 
supported on a finely grained silicon carbide plate, Saint- 
Gobain Ceramics, and covered by a platinum lid, the en- 
tire system is placed on a fire-brick. Fire-brick, platinum 
lid, silicon carbide plate and samples are cleaned with the 
same protocol : i) one ultra-sonic cleaning of fifteen min- 
utes in a millipore water/detergent (aquanox) 10% mix- 
ture, ii) two ultra-sonic rinsing of fifteen minutes in milli- 
pore water, iii) drying with dry nitrogen. After this clean- 
ing, all pieces are placed in a furnace at 60° C to dry the 
porous fire-brick. 



Annealing treatments on silica 

All silica samples were first submitted to a high temper- 
ature treatment: 1770 K during 30 minutes. This step is 
performed well above the Tq ~ 1450K to ensure full re- 
laxation of surface defects due to the polishing step. The 
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sample is then immediately taken out from the furnace; 
the platinum cap is only removed after 10 minutes in order 
to avoid dust in-layering on the molten surface 

Beyond this surface relaxation step, the silica samples 
have undergone different thermal treatments to test the 
dependence of surface roughness on thermal history. We 
followed annealing protocols proposed by Le Pare et al. 
|T3) who studied the effect of thermal history on the Ra- 
man spectrum of silica. The thermal treatments are listed 
in tabled Annealing times have been calculated to allow 
silica relaxation. After annealing each sample is removed 
from the furnace in the same condition as described above. 



Sample 


1500 


1350 


1200 


1100 


A 


30 min 








B 


30 min 


60 min 






C 


30 min 




120 min 




D 


30 min 






68 hours 



Table 1. Thermal history of samples 



Note that the long annealing time applied to sample 
D has caused a slight nucleation of a crystalline phase 
of silica under the surface. This is not expected to affect 
further interpretations as the crystals were subsurface. 

Post-annealing cleaning 

In order to eliminate contaminations due to storage and 
transport, samples were cleaned before AFM measure- 
ments. We checked that the surface was not altered by 
successive cleaning operations. The previous cleaning pro- 
tocol was used. A last step is added : a 45 minutes UV/O3 
cleaning to remove nanometric dust and also to obtain 
a hydrophilic surface improving surface imaging. Despite 
theses precautions; some dust was found on the surfaces . 
However, AFM scans show large dust free parts (typically 
15^m*15/im). 

Quantitative use of AFM for surface morphology. 

The surface morphology of all samples have been char- 
acterized by AFM height measurement in tapping mode 
(TM-AFM), using a Nanoscope III A from Digital Instru- 
ments with Al coated tip (BudgetSensors - model BS-Tap 
300 Al). The Al coating thickness is 30nm, the resonant 
frequency is 300 kHz and the stiffness constant is 40N/m. 
The tip radius is lower than 10 nm. The images have been 
recorded at a scan frequency between 0.8 and 1Hz for a 
resolution of 512 x 512 pixels. 

As our study is focused on the analysis of extremely 
smooth surfaces, special attention has to be paid to the 
data treatment. The roughness induced by capillary waves 
lies indeed in the sub-nanometer range and corresponds to 
the ultimate thermodynamic fluctuations of the interface. 
At this low roughness level, limitations of AFM can be 
observed in the bare data. 




I I : I 

mi^^E^^ 4096 nm 

-3 nm 4 nm 

1-a - Best paraboloid correction 





f M 4096 nm 

-1 nm 5 nm 

1-b - Correction by the best set of parabolas 
Fig. 1. An image of an annealed silica surface (A) acquired 
in the open air. a) correction by the best fitted paraboloid 
(p'rms = 0.43 nm) to compensate the systematic curvature ef- 
fect due to the spherical motion of the piezo-electric stage, note 
the residual bands due to a drift error of the tip; b) correction 
by the best set of parabolas (a rms — 0.25 nm). 



The sample is moved by a XYZ-piezo-electric stage 
and long profiles exhibit a slight spherical character. This 
systematic curvature in the height data can be eliminated 
by subtracting the best paraboloid to the surface. Fig. QJ 
-a) shows the result of this first operation. A remaining 
strong band pattern is observed; this results from a drift 
error of re-positioning of the piezo-electric stage. Strictly 
speaking AFM measurements thus have to regarded as a 
set of profiles more than an image and bare data must 
be corrected line by line. Fig.[IJ-b) shows a representative 
AFM image obtained after subtraction of the best set of 
parabolas (in a least mean square sense) to the set of pro- 
files. The statistical analysis is thus performed on these 
profiles. 
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4 Results 
4.1 Silica 

A set of AFM roughness measurements on the samples 
A,B,D has been made in order to determine the scal- 
ing behaviour of the roughness from the nanometer range 
up to the micrometer range. The typical experiment is 
the following : a first image is acquired over 8/xmx8/im, 
then a hierarchy of images of decreasing size (4/zm, 2/xm... 
down to 512 or 128 nm depending on the sample) is taken 
within this first image. We performed the measurements 
of 2 such series down to 128nm for sample A; 2 series 
down to 512nm for sample B and 3 series down to 512nm 
for sample D. On sample C we measured 14 images of 
2/imx2/xm and 7 sets of 2 imbricated images of respective 
size 8/imx8/im and 2/imx2/jm. 
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Fig. 2. Height variogram obtained from AFM measurements 
silica surfaces after various annealing treatments in the glass 
transition regime (symbols). The lines correspond to the pre- 
dicted logarithmic behaviour induced by the presence of frozen 
capillary waves (fit performed in the interval [16nm-500nm]) 
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Fig. 3. Height variogram obtained from AFM measurements 
on the two faces of a float glass, respectively in contact with 
liquid tin and with atmosphere (symbols). The fitted lines cor- 
respond to the logarithmic behaviour predicted by the presence 
of frozen capillary waves (fit performed in the interval [8nm- 
2000nm]). 



The statistical treatment consisted in computing the 
height variogram on each individual image and to average 
these data for the different images of a series. We thus 
obtain on Fig. [3 the evolution of the average quadratic 
height difference H(r) = (\h(x + r) — h(x)\ 2 ) x between 
points on the surface separated by a distance r. In the 
semi-log scale used in Fig. [21 we obtain a linear behaviour 
over more than 2 decades from a few nanometers up to a 
few micrometers. We thus obtain a very good agreement 
with the predicted capillary wave scaling proposed in Eq. 
(0. The values of the slopes obtained by a simple linear 
fit in the scaling region are reported in TableEl A physical 
discussion on these extracted parameters and the quanti- 
tative agreement with the theory is given in section [5] 



4.2 Float glass 

Using the same procedure as above, 6 sets of images of 
decreasing size from 8/imx8/xm down to 500nmx500nm 
were performed on the atmosphere and tin faces of the 
float glass samples. The atmosphere faces were more af- 
fected by residual dust pollution and only 3 of them could 
be quantitatively analysed. The height variograms are dis- 
played in semi- log scale in Fig. [3 Again we observe a loga- 
rithmic behaviour over more than two decades. The values 
of the slopes obtained by a simple linear fit in the scaling 
region are reported in Table [2J The tin face appears to be 
much smoother than the atmosphere face, which is con- 
sistent with the smaller interface tension expected for the 
glass atmosphere interface compared with the tin/glass 
interface. 



4.3 Statistical averaging, uncertainty 

To evaluate the accuracy of our results, we investigated the 
dispersion of the roughness measurements. This was done 
on the silica sample C. As a first test, 14 images of size 
2 //m were recorded within a primary zone of size 8 /xm. 
The mean value of the height variances obtained on these 
images was a 2 = 0.029nm 2 and the standard deviation 
0.005 nm 2 , that is to say 15 % of the mean value. As a 
second test, we performed 7 sets of 2 imbricated images 
of respective sizes 8 /xm and 2/xm measured at different 
places over the surface. For the 8 /xm scans, we obtain a 
mean value of 0.053 nm 2 and a root mean square of 0.007 
nm 2 that is to say 15 % of the mean value. For the 2 
/im we found a mean value of 0.030 nm 2 and a root mean 
square of 12 % of the mean value. We can conclude that we 
have roughly a dispersion of 15 % of the mean value. This 
gives us an estimate of the uncertainty on our roughness 
measurements and on the numerical parameters used to 
fit the logarithmic behaviour predicted by Eq. |5|. 



5 Discussion 

As shown in Fig[2] and [3 the height variograms show a 
clear logarithmic behaviour from the nanometer range up 
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Sample 


Tf(K) 


7 (J/m") 


a th (nm z ) 


a (nm 2 ) 


silica A 


1773 


0.3 


0.026 


0.022±0.003 


silica B 


1623 


0.3 


0.024 


0.020±0.003 


silica C 


1473 


0.3 


0.022 


0.020±0.003 


silica D 


1373 


0.3 


0.020 


0.019±0.003 


float atm. 


873 


0.35 


0.011 


0.024±0.004 


float tin 


873 


0.5 


0.008 


0.011±0.002 



Table 2. Estimation of the capillary wave slope parameter 
obtained by fitting the evolution of the mean square : H(r) — 
ar+b and comparison with reference values a t n = ^p- obtained 
from the annealing temperatures and estimates of the interface 
tensions. 



to the micrometric range. Table |21 shows a summary of 
the parameters characterising the frozen capillary wave 
regime. As references, we indicated the values kTp /tt^ cor- 
responding to the annealing temperatures used for silica 
and the temperature of the soda-lime glass at the end of 
the float tank; we used generally accepted estimates of the 
interface tensions silica/air, glass/air and tin/glass at high 



Several comments can be made upon these results. 
First, within experimental errors we obtain a clear log- 
arithmic scaling of the spatial correlation with physical 
parameters consistent with the scenario of capillary waves 
freezing. 

In this naive scenario we assume that the surface freezes 
at a well-defined temperature Tp. Some more complex 
phenomena may characterise the freezing step, in particu- 
lar one can think of a dependence of the freezing temper- 
ature on the spatial frequency dependence of the capillary 
waves. However such glassy behaviours, leading to depar- 
tures from the simple logarithmic scaling could not be 
identified within experimental errors. 

The logarithmic behaviour could be observed over 2 
to 3 decades from the nanometer range to the micrometer 
range. Our data do not show any clear upper cut-off. Our 
AFM set-up being limited to a maximum scan length of 
17 /urn, larger scales were not explored. Note however that 
the pollution of the surfaces by dust or impurities makes 
it more and more difficult to analyse the data beyond 10 
/nm. As discussed above, the lower cut-off of the logarith- 
mic range is supposed to correspond to an atomic length 
scale. In our case, this cut-off has an instrumental origin 
and amounts to a few nanometers, which correponds to the 
AFM tip resolution. Note again that the logarithmic scal- 
ing induces a divergence for both small and large scales. 
Variations of the lower cut-off value due to a change of 
rounding of the AFM tip thus may induce a vertical shift 
of the roughness data (see Eq. J5||). 

In this framework of frozen capillary waves, the slope 
measured in the semi-log plot is directly related to the 
ratio Tp/') where Tp is the freezing temperature and 7 
the interface tension. 

In the case of silica we observe first a very good agree- 
ment between experimental slopes and references values. 
We observe a very low growing tendency of this slope with 
the annealing temperature. This may thus correspond to 



a slight increase of the effective freezing temperature but 
this tendency is certainly not clear enough to give any 
firm conclusion. The typical slope obtained around 0.20 
nm 2 are very close to the value 0.22 nm 2 corresponding 
to the glass transition temperature of silica Tq — 1450K 
and an interface tension 7 = 0.3J.m~ 2 . Let us note in 
addition that, because of thermal inertia, it is not surpris- 
ing that for the annealing processes at high temperature 
the effective freezing temperature (the experimental slope) 
be lower than the annealing temperature (the reference 
slope). 

In the case of float glass we observe a striking con- 
trast between the two faces. The tin face is noticeably 
smoother than the atmosphere face as can be expected 
from the contrast of interface tensions. Focusing on the 
quantitative aspect, we note that the experimental slopes 
are quite higher than the reference values, particularly in 
the case of the atmosphere face. This may obviously be due 
to an underestimation of the freezing temperature and/or 
an overestimation of the interface tensions. Let us insist 
on the latter point, the experimental measurement of in- 
terface tensions at high temperature (sessile drop methods 
and further evolutions) remains delicate and large varia- 
tions of the measured values can be found in literature|15j. 
The interface tension can indeed be very dependent on the 
atmosphere and on the composition. 



6 Conclusion 



Performing AFM roughness measurements on annealed 
silica surfaces and on float glass surfaces we could con- 
firm that the roughness of "fire polished" glasses can be 
quantitatively described by frozen capillary waves in the 
micrometric range. In particular we could identify a clear 
dependence of the roughness spectrum on the interface 
tension and we obtain indications of a slight dependence 
upon the freezing temperature. At length scales below the 
capillary length, the industrial float process thus allows 
to produce flat glass surfaces characterised by the lowest 
thcrmodynamically possible roughness; the only way of re- 
ducing these fluctuations would consist of increasing the 
interface tension or of decreasing the "freezing" temper- 
ature of capillary waves. In addition, we could show that 
the description of glass surfaces by capillary waves gives 
access to the ratio freezing temperature on interface ten- 
sion Tp /7 and can thus be considered as a complementary 
indirect way of measuring the two quantities. 
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